Separation of time-scales and reparametrization invariance for aging systems 
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We show that the generating functional describing the slow dynamics of spin glass systems is 
invariant under reparametrizations of the time. This result is general and applies for both infinite 
and short range models. It follows simply from the assumption that a separation between short 
time-scales and long time-scales exists in the system, and the constraints of causality and unitarity. 
Global time reparametrization invariance suggests that the low action excitations in a spin glass 
may be smoothly spatially varying time reparametrizations. These Goldstone modes may provide 
the basis for an analytic dynamical theory of short range spin glasses. 
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Much progress has been made in recent years in under- 
standing the non-equilibrium dynamics of glassy systems. 
Late after a quench to low temperatures, two distinct 
dynamic regimes develop [1,2]. At short time-differences 
the dynamics is similar to that in the disordered phase. 
The correlations depend only on time differences [time 
translation invariance (tti)] and decay towards a non- 
vanishing Edwards- Anderson (ea) order parameter. The 
correlations and their associated responses are related by 
the fluctuation-dissipation theorem (fdt). At long time- 
differences, for a freely relaxing system, physical quan- 
tities relax very slowly and depend on the waiting time, 
so that correlations depend on two times rather than on 
time differences. The separation of time-scales exists also 
when a glassy system is gently driven, whether the full 
dynamics becomes stationary or not. In the slow regime, 
the fdt is replaced by an out-of-equilibrium fluctuation 
dissipation relation (oefdr) between each bulk response 
R and its associated bulk correlation C: 

R(h,t 2 ) = 0™AC(h,t 2 )} J-C(ti,t 2 ), (1) 

Ot 2 

where ti > t 2 , and /3 EFF [C] is the inverse effective temper- 
ature measuring deviations from the fdt (e.g. [3]). These 
properties, observed in experiments and simulations [1], 
are embodied in the concepts of weak ergodicity break- 
ing and weak long term memory [2,4,5]. This scenario 
has been demonstrated analytically for mean-field glassy 
models, such as the Sherrington-Kirkpatrick model [5]. 
However, for short-range models such as the Edwards- 
Anderson model [6], no analytic solution is known for 
either the statics or the dynamics. The work here pre- 
sented uncovers a symmetry that strongly constraints the 
properties of such a dynamic solution. 

The oefdr is consistent with the observation that the 
equations of motion of a large class of glassy models are 
invariant under time reparametrizations, t — > h(t) in the 
slow regime [3,5,7-9]. In particular, 



R(ti,h 



dh 



R(h(h), h(t 2 )), 
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are related by the same oefdr as in Eq. (1). Equa- 
tions (2) and (3) are particular cases of the general cor- 
relator Reparametrization Group (RpG) transformation 



dh 



dh 
di 2 



G(h(h), h(t 2 )) , (4) 



with A§ and defined in Ref. [8] as, respectively, the 
advanced and retarded scaling dimensions of G under 
the reparametrization t — > h(t) of t\ i2 . In particular, 
, = 0, and Ag = 1 [see Eqs. (2) and (3)]. 
In this paper we show that: i) RpG invariance exists 
at the level of the generating functional for the long- 
time dynamics, and ii) it holds for short-range models 
and their local, and not only bulk, quantities. These are 
much stronger results than those of Refs. [3,5,7-9], which 
only applied to mean field models and were restricted to 
global quantities. Furthermore, the existence of a con- 
tinuous symmetry in the action allows us to predict the 
presence of a Goldstone mode controlling the dynamics, 
which may be the basis for developing a systematic an- 
alytical theory of dynamical fluctuations in short range 
spin glasses. 

In order to prove our result, we take three steps: i) we 
determine the long-time action by using the Renormal- 
ization Group (RG) in the time variables, assuming that 
there is a separation between short and long time scales, 
ii) we analyze the surviving terms in the action, showing 
that they are RpG invariant, and iii) we show that the 
measure in the functional integral is also RpG invariant. 

For the sake of concreteness we discuss the general dis- 
ordered spin model defined by 



(5) 



C{t 1 ,t 2 )=C(h{t 1 ), h{t 2 )) 



(2) 



where H\ oc contains arbitrary self-interactions that, for 
example, place soft or hard (±S) constraints on the 
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spins. Typically, the random couplings are Gaus- 
sian distributed according to P(Jij) oc exp[—Jfj/(2Kij)], 
with Kij the connectivity matrix. In the EA model, 
Kij = J 2 1 z if i,j are nearest neighbors or zero other- 
wise (z is the coordination of the lattice). We consider 
a quantum extension [8,10], with quantization rules im- 
posed via a kinetic term, ^2f = i pf / (2M) , where the mo- 
menta pi are canonically conjugate to the "coordinates" 
Si, \pi,Sj] — —ihSij. The use of bosonic variables to rep- 
resent the spins (e.g. quantum rotors, which lack Berry 
phases) simplifies the presentation considerably although 
the arguments below can be extended to other glassy sys- 
tems such as the SU(N) Heisenberg model [11]. The sys- 
tem is linearly coupled to an equilibrated environment. 

We study the dynamics with the Schwinger-Keldysh 
closed time-path formalism [8,10-14] in which the spin 
variables acquire another index labeling the two Keldysh 
branches, Si — > Sf, a = 0,1. We work in the rotated 
basis a? = &i = (S? - Sf)/V2 and a\ = a % = (Sf + 
Sf )/V2. This notation renders the classical limit, % — > 0, 
more transparent [10]. 

Once the disorder has been integrated out, we intro- 
duce a local Hubbard-Stratonovich field Qi(t\, t 2 ) [15] 
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that decouples the four spin term generated by the av- 
erage over the Jy. The choice of indices 0, 1 anticipates 
the result that we prove, namely invariance of the action 
under the transformation Qf b (ti,t 2 ) — ► Qf b (ti,t 2 ), with 



(dh_ 



dh_ 
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Qf(ft(*i),M*2)) (7) 



[cf. Eq. (4)]. This Ansatz for how the Q fields transform 
is motivated as follows. The Q K component is a time- 
dependent local measure of freezing (analogous to the 
Edwards- Anderson order parameter for the static case) 
since it is related to the correlator at the same site but 
at two different times; hence the choice of vanishing di- 
for this slowly decaying quantity. 
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mensions A 

The scaling dimensions for the retarded and advanced 
components Q R A are suggested by the fdt, once the 
dimensions of Q K are fixed. 

In the usual Keldysh approach, one keeps only the 
three non- vanishing components, K, R, A, since Keldysh 
propagators are related to expectation values (Qf~), (Qf) 
and (Qf), while (Qf) = [12,13]. Here we include the 
fourth component because we are also interested in the 
fluctuations of the Qi's. The generating functional reads 

Z = J [DQ] cxp (-S K [Q\ - Snl[<2]) , (8) 

S K [Q] = \J2 k ^ [ QtihM) Qf(t 2 ,h), 



S Nh [Q] = - m f[Do-}[Da] cxp {iS[a, a; Q] + iS SPm ) , (9) 
S[a,a;Q] = J2 J dtldh <^ Qffah) ■ ( 10 ) 

i,a,6 

^spin includes all the local spin dynamics, including the 
kinetic term, the self-interactions in H\ oc , and the cou- 
pling to the bath. The overline is a shorthand notation 
such that = 1 and 1 = 0. All integrals start at t = 0. In 
order to observe noncquilibrium glassy dynamics, we first 
take the thermodynamic limit, N — > oo, consider times 
such that the upper limit in the time-integrals diverges 
subsequently [3], and use a weak coupling to the bath. 
Upon integration over the spin variables erf in Eq. (9): 

oo . 

Snl [Q] = X) -, E dtl --- dt2n G -V;:::X (*i - ■ • • 

n—l ^l : ...,i n 

x QtrihM)-- - C"" 2 "^-!.^) - (11) 
with Gi\ 7 - : x n (h,...,t 2n ) «(ti)---(C»(t2n)>. 

[( • ) denotes an average over the fields a, u with the 
weight exp(iS SP i N )]. The non-linear terms in Eq. (11) arc 
generated because H\ oc is non-quadratic. This expansion 
contains only connected terms in the sense that none of 
the powers of Q can be factored. Each Q serves as a 
source for a pair of cr's, hence the 2n-point correlations 
q(2u) arc not ncccssar iiy connected in terms of individ- 
ual g 's. In general, the contribution at any order n in Q 
can be evaluated by rewriting the G^ 2 "^ into sums over 
products n^Li Gc 2m °' of connected Green's functions 
Gc 2m "\ with mi + • • • + rtiM =n. An example, to order 
Q 3 , is shown in Fig. 1; notice that the three Q's form a 
connected bubble diagram, though at the same time it 
contains both a 2-point and a 4-point connected Green's 
function for the cr's. Later on we will use this example 
to illustrate some of the steps in our arguments. 



X:- 1 ' -a, (t 4 

a-| \ >2 

o i (t,) 

CT u (t e )y. 3 

3 



QAh.U 



Qi (t 3 ,t 4 ) 



FIG. 1. A cubic term contributing to Snl[Q]. The dashed 
lines are a's, which together with the shaded regions represent 
connected Green's functions for the cr's. 

i) RG in the time variable - We assume there is an 
initial short-time cutoff To = f2 _1 . We then separate the 
fast and slow components of the field Qi, 
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Qf(ti,t 2 ) 



Qf FAST (h,t 2 ) if r < |*i - t 2 \ < bT , 

SLOW 



(ti,t 2 ) if br < \h -t 2 \ 
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for a rescaling b > 1 (and 51 = In b > 0). The 
couplings that flow in the RG are the coefficients 

o 2 n(* 1 ' • • • ' * 2 «) * n -^l- (H)- ^ s USU& 1; eacn 



transformation involves two operations. The first is the 
integration over the fast modes <3i FAST , which we per- 
form by representing S NL as a path integral over spin 
variables erf [Eqs. (9) and (10)]. The second is rescaling 
time t — > t/b to restore the cut-off to its original scale 
r = fi -1 , accompanied by a rescaling of the Q fields, 
performed on the expression in Eq. (11). Each of these 
operations generates a change of order 51 in the action. 
The quadratic term in the action, Sk[Q], plays a role 
analogous to the kinetic energy term in a usual RG cal- 
culation: it does not mix fast and slow modes, and the 
rescaling of fields is chosen so as to keep it invariant. 

The integration over fast modes yields a non-local four 
spin interaction similar to the one obtained by perform- 
ing the disorder average, but fundamentally different in 
that it is short ranged in time. This extra spin-spin inter- 
action leads to a change in the 2n-point spin correlation 
functions, 5G {2n ^ = 51 [dG^/ W]| fast . 

As mentioned above, each coefficient G^ 2 ") can be ex- 
pressed as a sum of products of connected Green's func- 
tions. Consider one of these connected Green's functions, 
which plays the role of a coupling for the Q's. In the same 
way that in a usual RG calculation (in real space) nonlo- 
cal interactions flow into local interactions (i.e. the only 
couplings that are left are those connecting fields at the 
same point) , in our case all couplings connecting different 
times flow into couplings that only connect equal times 
(i.e. other terms are irrelevant compared to these terms). 
Using this property, the term in Fig. 1 reads: 

fdr G^% (r, 0) J d?r G C ££™ B (V, r", r'", 0) 

x J dhdh Ql ia2 (h,t 2 ) Ql 3ai (t 2 ,t 2 ) Qg-fah) . 

Here and in the next equation, the r variables represent 
time differences (time can be shifted because of tti) and 
their range of integration is extended to be (— oo,+oo). 
For a general term in the slow action, we obtain an in- 
tegral of a product of Q fields over time variables (one 
time variable per connected spin Green's function), with 
prefactors given by the integrals of the connected spin 
Green's function over all possible time differences that 
define the coupling constants: 



H,...,»2m 
JOl,...,02ti 



J d 2m -\ 



G 



cai,...,a 2 , 



(n, . . . , r 2m _i, 0) 



= iim xcaV;::.',a m 2m (wi, ■ ■ ■ , w 2m -i) . (12) 

They correspond to physical zero-frequency (dc) gener- 
alized correlators \c C °f the spin variables. 

Let us now turn to the effect of restoring the cut-off to 
the original scale by rescaling the times, t —* t' = t/b = 



te- 51 , and the fields, Q* ia2 {t u t 2 ) -» Q'T^iAA) = 
b a 1+ a 2 Q»i»2(t 1)f2 ). it i s easy to see that it satisfies 
two conditions: it produces the same effect as an RpG 
transformation with h(i) = bt [cf. Eq. (7)], and it leaves 
the quadratic term Sk invariant. In the example: 

/«x[Q] = J dhdt 2 QT^itiM) Qi 3ai (t2,h) QT 3 a6 {t2M) 
= b 2 J dt[dt' 2 b- < - ai+a6 H- < - a2+a3+ai+a ^ 

xQt a2 (*i>4)Q'r 4 (4,4)QT 6 (4,ii) 

_ ^l-(ai+o 6 ) ^l-(a2+a 3 +a 4 +a 5 ) j \Q'] 



In a general term, for each time variable t a there is a 
coupling g a = ga~i',"!*a£? a , and also an exponent A Q = 
ai + - ■ -+a 2ma > originating from the rescaling of fields. 
In the example A x = a\ + a e and A 2 = a 2 + a 3 + a A + a 5 . 

Combining the contributions from both the integration 
of fast modes and the rescaling, we obtain the following 
flow equation for the coupling constant g a : 



dg a /, a n , d 9a 

^ = (l-A a ), Q + — 



fast 



(13) 



Terms with A a = occur only when all the a's are 0. 
These terms would be naively relevant (according to their 
engineering dimensions), but in fact they vanish identi- 
cally due to the constraints of normalization and causal- 
ity for G { c n) [12,13]: 



/-Y»lr">»2n a I, 

^0,0, ...,0 I 1 !' 



i hn a ) = (*l) ' 



! (*2n a )> = 0. 



There remain only terms which are either marginal (for 
A a = 1) or irrelevant (for A a > 2) according to their en- 
gineering dimensions. Under the assumption of a sepa- 
ration between short and long time scales, at some point 
in the RG flow all of the fluctuations associated with 
short time scales will have been integrated over, (but the 
flow will not have yet reached a point were the long time 
scales are probed) and therefore any new integrations of 
fast modes produce no change in the coupling constants, 
i.e. the second term in the r.h.s. of Eq. (13) is zero. This 
assumption is used in an analogous way in the solution 
of mean field spin glass models [1-3]. Physically, the sep- 
aration between short and long time scales correspond 
to a saturation of the generalized DC correlators for the 
spins Xc° defined in Eq. (12) to a finite value. Therefore, 
the existence of a separation of time scales implies that 
the engineering dimension actually determines the long 
time behavior, and only marginal (A a = 1) terms are left 
in the effective long-time action. 

ii) RpG invariance of the effective long-time action. 
For a RpG transformation applied to a generic term 
in 5 NL , each integral over time gives rise to a factor 
(dh(t a )/dt a ) a , and since only terms with A a = 1 are 
present, the derivative factor yields exactly the Jacobian 
necessary to make the term RpG invariant: 
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'T\dt 



dh 
dt a 



-/n*.(^)--/n*.-- 

To illustrate this point, let us perform a RpG transfor- 
mation on our example term: 



IM = / dhdt 2 Q% a2 (ti,h) Q? 2 a *{t 2 ,t 2 ) QT'^M) 



dhdt 2 (9/i(ti)/ati) ai+a6 (dh(t 2 )/dt 2 ) a2+a3+a4+a5 



-S 

x Q^(h(h),h(t 2 )) Q^(h(t 2 ),h(t 2 )) QTf^Kh)^^)) 
= J dh x dh 2 Q^ a2 (h u h 2 ) Q^ ai (h 2 ,h 2 ) Q? 3 5 ° 6 (/i 2 ,/u) 

- ; 

where we have used that a\ + a 6 = Ai = 1 and a 2 + a 3 + 
d4 + as = A 2 = 1. For the Sk term, the same argument 
applies with A = 1 replaced by a + a = 1 and b + b = 1 . 
Therefore the long-time action is RpG invariant. 

iii) RpG invariance of the measure - Under the RpG 
transformation of Eq. (7), the Jacobian for the functional 
integral over the Q fields is simply: 



J 



DQ 
DQ 



nn 

x ti,t 2 



dh dh 
dhdh 



J d d xd tl dt 2 in | 



The Jacobian depends on h(i), but not on the fields Q, 
and therefore the generating functional is RpG invariant. 

The RpG invariance implies that the action describing 
the long time slow dynamics of a spin-glass is basically 
a "geometric" random surface theory [16], with the Q' s 
themselves as the natural coordinates. The original two 
times parametrize the surface. Physical quantities, as the 
bulk integrated response x(ii>*2) = J t 1 dt' (Q R (f ,t 2 )) 
and correlation (Q K )(ti,t 2 ) have scaling dimension zero 
under t — > h(t) [5,8] as well as their local counterparts 
that are directly related to the Q^s. A possibly related 
gauge-like symmetry has also been noted in the replica 
approach [17]. 

The emergence of the RpG invariance may provide a 
novel, completely dynamical, angle to address the still 
poorly understood short-range spin-glass problem ana- 
lytically. We have shown that the global reparametriza- 
tion, t — > h(t), is a symmetry of the slow dynamical ac- 
tion. The particular scaling function h(t) selected by 
the system is determined by matching the fast and the 
slow dynamics. It depends on several details - the exis- 
tence of external forcing, the nature of the microscopic 
interactions, etc. In other words, the fast modes which 
are absent in the slow dynamics act as symmetry break- 
ing fields for the slow modes. The global RpG invari- 
ance of the slow action suggests that the low energy 
physics of the glassy phase could be described by slowly 



spatially varying reparametrizations t — > h(x,t). Basi- 
cally, we propose that there are Goldstone modes for the 
glassy action which can be written as slowly varying, spa- 
tially inhomogencous time reparametrizations. Compar- 
ing with the O(N) non-linear sigma model [18], global 
time-reparametrizations are analogous to uniform spin 
rotations, while local t — ► h(x, t) reparametrizations de- 
scribe the spin waves (fluctuations on the uniform solu- 
tion) . Numerical tests in the 3D EA model are consistent 
with this conjecture [19]. 
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